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ABSTRACT 

Baryonic acoustic oscillations (BAOs) modulate the density ratio of baryons to dark 
matter across large regions of the Universe. We show that the associated variation 
in the mass-to-light ratio of galaxies should generate an oscillatory, scale-dependent 
bias of galaxies relative to the underlying distribution of dark matter. A measurement 
of this effect would calibrate the dependence of the characteristic mass-to-light ratio 
of galaxies on the baryon mass fraction in their large scale environment. This bias, 
though, is unlikely to significantly affect measurements of BAO peak positions. 
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1 INTRODUCTION 

The rapid acoustic waves in the radiation-baryon fluid prior 
to cosmological recombination were not followed by the dark 
matter at that time. Following recombination, the baryons 
were freed from the strong radiation pressure and fell into 
the gravitational potential fluctuations of the dark matter. 
As a result, the fractional difference between the density 
fluctuations of baryons and dark matter decreased steadily 
with cosmic time. But since the baryons amount to a sizable 
fraction of the total mass density of matter (f2(, /Q m ~ 17%), 
the gravitational effect of the baryons on the dark matter im- 
printed baryonic acoustic oscillations (BAOs) on the matter 
power spectrum. The characteristic comoving scale of BAOs 
~ 100 Mpc (corresponding to the sound horizon at recombi- 
nation), provides a yardstick that can be used to measure the 
dependence of both the angular diameter distance and Hu b- 
ble parameter on redshift (see review bv lEisenste"mll2005l ). 

When analyzing galaxy surveys, it is often assumed 
that galaxies ar e biased tra cers of the underlying mat- 
ter distribution (|Kaiser| [l984). with a bias factor that ap- 
proaches a constant value on sufficientl y large scales where 



density fluctuations are s t ill linear (e.g..lMq fc Whitdfl996l ; 
iTegmark fc Peebles! 1 19981 ; ISheth et afll200ll) . However, the 
imprint of primordial acoustic waves on the baryon fluid 
at recombination introduced a scale-dependent modulation 
of the ratio between the density fluctuations of baryons 
and dark matter that has not been completely erased by 
the present time. A large-scale region with a higher baryon 
mass fraction than average (in the perturbations that lead 
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to galactic halos) is expected to produce more stars per unit 
total mass and hence result in galaxies with a lower mass- 
to-light ratio. 

In this paper we characterize the associated scale- 
dependent bias in flux-limited surveys of galaxies. The ratio 
between the power spectra for fluctuations in the luminosity 
density and number density of galaxies is expected to show 
BAO oscillations that reflect the large-scale variations in the 
baryon-to-matter ratio. 

In §2, we formulate the oscillatory BAO signature on 
galaxy bias in terms of a simple analytical model. The quan- 
titative results from this model are presented in §3. Finally, 
we summarize our main conclusions in §4. 



2 THE MODEL 
2.1 Basic Setup 

Since galaxies sample the high peaks of the underlying mat- 
ter density, they are biased tracers of the matter density. 
When the clustering of galaxies is usually analyzed, the 
bias is considered simply with respect to the matter density, 
without separating out the effects of the baryons. As long 
as the baryon fluctuations follow the same spatial pattern 
as that of the dark matter, biasing with respect to each of 
them cannot be separated since this separation is degenerate 
with an overall change of the bias factor, which is not known 
apriori. However, since the BAOs induce a scale-dependent 
difference between the baryons and dark matter, it becomes 
important to consider their influence on galaxies separately. 

Consider the power spectrum of fluctuations in the 
galaxy number density n ga i and in the luminosity density 
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Pl- For a given galaxy population, 

Pl = n gal x (L) 



(1) 



where (L) is the mean luminosity of the galaxies. Since 
galaxy formation is driven by halo collapse, which depends 
on the evolution of the overall matter perturbations, the 
number density fluctuations <5 n are driven by the fluctuation 
5tot in the total matter density, with a bias b n that should 
be approximately constant on large scales (for a fixed galaxy 
population) : 



S n = b n 5t, 



(2) 



The mean luminosity of galaxies may depend on their envi- 
ronment through their merger rate history, which is corre- 
lated with the local matter density. This can lead to fluctu- 
ations <5l in pl with a different bias factor that should also 
approach a constant on large scales: 



(&n + &L;t)<5tot 



(3) 



where the overall bias factor of the luminosity density with 
respect to the total matter includes the number density bias 
b n as well as a possible additional bias &L;t from the depen- 
dence of (L) on the matter density. 

However, the luminosity is also affected separately by 
the baryon fluctuations, since the luminosity depends on 
the gas fraction in halos /b. Regions that have halos with 
a higher baryon fraction will proportionally have more 
baryons in the galaxies within them. If, e.g., we assume that 
the star formation rate per baryon is on average constant, 
then (L) oc /b. In fact, the dependence of the luminosity 
on the gas fraction is likely to be non-linear. For instance, 
in sim ple models for disk formation within halos l|Mo et al.l 
1998), the disk radius is approximately independent of the 
gas fraction. Thus, if we assume that the disk mass is a fixed 
fraction of the halo gas mass, then the typical gas surface 
density within the disk varies in proportion to the overall 
halo gas fraction. Acc ording to the Schmidt-Kennicutt law 
(e.g.. lKennicuttlll998l ). the star formation rate in the disk 
should vary with the gas surface density to the power 1.4. 
Thus, in general we assume that 



(L) «(/ b ) 6 



(4) 



where these simple considerations suggest that i>L ; f ~ 1-4. 
The notation for this power index is chosen since equa- 
tion @ (together with equation [3} implies fluctuations 



5l = (&n + &L;t)(5tot + &L;f5f , 



(5) 



where <5f is the perturbation in the halo gas fraction /b. 
Thus, &L;f is the bias factor of the luminosity density with 
respect to the halo baryon fraction. Note that in our notation 
all the perturbations are the actual ones at the considered 
redshift (i.e., we do not use the common practice of linear 
extrapolation to redshift zero). 

2.2 Halo Baryon Fraction 

We would expect the baryon fraction within halos to reflect 
that of their surroundings, but the precise relation is com- 
plex due to the non-linear process of halo collapse. Here we 
employ reasonable simplifications to derive an approximate 
analytical result, which is partly verified and quantified by 
simulation results shown in § 3. 



We find it useful to analyze the baryon fraction in sev- 
eral steps, where the first step is to avoid halo collapse and 
simply consider 



7b 



Pb 

ptot 



(6) 



where we use 7 b for the general baryon fraction and reserve 
fb for the baryon fraction inside halos. The mean of this 
quantity is the cosmic mean baryon fraction: 



7b 



and its fluctuation is simply 



Ob 



rSt, 



(7) 



(8) 



Here we have measured the fractional difference between the 
baryonic and total matter fluctuations with r = (<5b/<5tot) — 1, 
in general a function of both wavenumber k and redshift. 

In reality, halos form out of perturbations that eventu- 
ally grow to an overdensity of hundreds, making the contri- 
bution of the mean density negligible, and thus we expect 
the baryon fraction to reflect the relative mass of the baryon 
perturbation that formed the halo: 



Q b S b - S b 
Jb = t^— ? — = 7b t — 

KraCtot Otot 



(9) 



Before discussing non-linear collapse, we wish to apply this 
equation to the linear perturbations that will form a halo, 
but even in the linear case we cannot easily apply this equa- 
tion in Fourier space, since halos form out of a sum of per- 
turbations on all scales, and taking a ratio as in equation @ 
is a non-linear operation. 

To make further progress, we mak e a separation of scale s 
(also called a peak-background split: [Cole fc Kaiserlll989h . 
where we assume that the fluctuations that we wish to ob- 
serve (in the galaxy luminosity, etc.) are on much larger 
scales than the (initial comoving) scales that formed the ha- 
los. Typically, we are interested in measuring fluctuations 
on BAO scales, which are ~ 2 orders of magnitude above 
the halo formation scale of galaxies. Thus, we separate out 
the linear halo perturbations (i.e., the initial perturbations 
that will form a halo, linearly extrapolated to the formation 
redshift of a particular halo): 



Stot — S tot + St ot , 
Sb = 5 l b + 5 b = (1 + n)$ 



+ (1 + r s )St 



(10) 
(11) 



where the relative difference between the baryonic and total 
matter perturbations is r\ and r B on large and small scales, 
respectively. 

We now use the standard result of spherical collapse, 
that a forming halo has a linear S to t = S c , where the critical 
density of collapse 8 C is independent of mass (and equals 1.69 
in the Einstein de-Sitter limit, valid over a wide range of red- 
shiftfS). We also assume that we are considering sufficiently 
large scales so that S l tot can be treated as a perturbation of 



1 The value of 8 C decreases at low redshift due to the cosmolog- 
ical constant, and at very high redshift due to the effects of the 
bar yons and radiation. How ever, at all z < 20 the change is below 
1% jNaoz fc Barkanall2007f) . 
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Stot (or <5t t), and that r\ and r s are also small quantities. 
Then the mean baryon fraction in halos is 

/b=7b(l + r s ), (12) 

and the lowest order perturbation is derived to be 



5{ 



n 



5c 



— Ot.n 



(13) 



We now use the actual value of r(k) (see § 3), specifi- 
cally the fact that it approaches a constant on scales be- 
low the BAOs, with a value (depending on redshift but not 
k) that we denote tlss (fo r Large Scale Structure) follow- 
ing |NaozTBarfen3 (|200Tt) . Thus, in the just-derived equa- 
tions we can treat r s = t-lss as a constant (at a given red- 
shift), since most of the density 5 C needed to form a halo 
comes from scales well below the BAO scale. Thus, the mean 
baryon fraction in halos is 



fb = 7b (1 + »"lss) i 
while on large scales (i.e., small k) the fluctuation is 
r(k) 



St = 



r L ss 



(14) 



(15) 



The remaining issue is the effect of non-linear collapse, 
and the relation between the baryon fraction in the linearly- 
extrapolated halo perturbation and the baryon fraction in 
the actual virialized halo. We show simulation results in § 3 
that only test the mean baryon fraction in halos (i.e., equa- 
tion I14[) but do so over a range of redshifts, and suggest 
that halo collapse enhances the effect and results in an ef- 
fective value of tlss that is amplified by a factor of several. 
One way to understand this enhancement is to consider the 
variation of 9"lss with time. It declines (in absolute value) 
approximately as r oc 1/a (where a = 1/(1 + z) is the scale 
factor), since (5tot —Sb) ~const while <$ t ot oc a (until the cos- 
mological constant becomes significant at low redshift). The 
decline of t-lss with time is of critical importance, since we 
are computing it according to linear theory, and it may not 
be appropriate to extrapolate tlss all the way to the halo 
formation time when we evaluate it in equation (|14[) . The 
baryon fluctuations, which were erased on small scales be- 
fore cosmic recombination, later continuously catch up with 
the dark matter (and thus with the total matter as well) 
in linear perturbation theory. However, once a perturbation 
begins to form a halo and enters the non-linear stage of col- 
lapse, we expect that the rapid collapse will bring with it 
only the baryons already present within the perturbation, 
and the continued decline of the linear-theory tlss will be- 
come irrelevant for the halo gas content. The upshot is that 
the simulations suggest that if we use the linear-theory n,ss 
(and, we assume, more generally for r(k)) then we must 
multiply it by an effective amplification factor A r : 



fb = 7b(l + A r 7-LSs) , 

Sf = 4 L ( r ( fc ) - r LSs]d"tot 

Oc 



(16) 
(17) 



The resulting fluctuations in the luminosity density (equa- 
tion [S| are 

Sl = (6n + &L;t)d"tot + &L;A[r(ft) - 7'LSs]o~tot ) (18) 



where 



&L;A = &L:f 

' Oc 



(19) 



is an effective bias factor that measures the overall depen- 
dence of galaxy luminosity on the underlying difference A 
between the baryon and total density fluctuations. 

2.3 Flux Limits 

We have assumed thus far that we observe a fixed galaxy 
population, regardless of the varying luminosity of its mem- 
bers. In reality, observed samples are limited by flux, or 
equivalently by luminosity if for simplicity we consider 
galaxies at a single redshift. Suppose the fraction of galaxies 
above luminosity L is 



F(L) = / 4>{L')dL' 

J L'=L 



(20) 



where <j) is the luminosity function. Then the observed num- 
ber density of galaxies is 



n bs = n ga iF(L) , 
and the luminosity density of these galaxies is 
Pl = n sa x(L)F(L) , 

where 

1 



(L) 



F(L) 



L'<t>{L')dL' 



(21) 



(22) 



(23) 



We assume for simplicity that the same luminosity dis- 
tribution holds in different regions, except that the lumi- 
nosity of all galaxies is enhanced or diminished uniformly 
in response to changes in the total density and the halo 
baryon fraction, as discussed in § 12.11 If a sample only in- 
cludes galaxies above a detection threshold Lmin, then we 
can analyze the variations of F{L) by keeping <f> fixed and 
varying the effective threshold L m i n , while in pl we also 
include the perturbation in the luminosity of each galaxy. 
From equation (120 ^ we obtain a relative fluctuation 

8p = OninlfeLitOtot + &L;f<5f] , (24) 

where the dimensionless coefficient 

0(Lmin) , _n 

6min _ F(L min ) ■ (25) 

The dependence of luminosity on the halo baryon frac- 
tion introduces a dependence of the galaxy number density 
on the baryon fluctuations (i.e., on r(k)). Putting our results 
together, for a flux-limited survey we find 



Sn — (fen + C m in&L;t)0~ti 

and 



+ C m infeL;A[r(fe) - r L ss]<5tot , (26) 



Sl = [fen + (l + Dmin)&L;t]Jtot + (l + -Dmin)bL;A[r(A;)-rLSs]5tot , 

(27) 

where 



(28) 



with (L) evaluated for L — L m i n . 

In the limit where L m i n is well below the peak of the 
luminosity function, C m i n and D m in both approach zero, 
and these expressions simplify to the previous ones (equa- 
tions [2] and [18]). In the opposite limit, e.g., in the exponen- 
tial tail of a Schechter function, we can approximately set 
4>(L) oc e~ L / L *, and then C m i n = L m i n /L, and D m i n = 
Cmi n Lmin/(Lmin + L») are both > 1 when L min > L*. 
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2.4 Observational Proposals 

As we have shown, both the galaxy luminosity density and 
(for a flux-limited sample) number density depend on the 
halo gas fraction. The scale-dependence of the relation be- 
tween the baryon and dark matter fluctuations implies that 
the BAOs can be observed in ratios that previously would 
have been expected to be scale-independent. 

One proposal is to compare the power spectrum of fluc- 
tuations in the galaxy number density (P n ) with that of the 
luminosity density (Pl), with both measured for the same 
galaxy sample. Taking the ratio may help to clear away some 
systematic effects that affect both power spectra. Their ratio 
(square-rooted) should have the form (assuming r(k) <IC 1): 



1/2 



-±L\ = B 1 {l + B 2 [r(k)-r hss ]} 



(29) 



where the various bias factors enter into the coefficients B\ 
and £>2- If we denote the bias ratio 6 r = &L ; t/&n, then 



Si = 



and 



B 2 = 



1 + (1 + J> min )br 
1 "h C'min^r 



l + D min Cn 



(30) 



(31) 



(1 + CmiA) • [1 + (1 + Anin)&r] ' 

Note that in the limit where most of the galaxy popula- 
tion is observed (i.e., the flux limits are unimportant), these 
expressions simplify to Bi = 1 + b T and B 2 = &l ; a/(&ii-Bi). 

In practice, using these expressions is not as daunting 
as it may appear. For a given galaxy sample, C m i n and -D m i n 
can be calculated from the measured luminosity function. 
This leaves two unknowns, b r and the ratio &L ; A/&n- Within 
the ratio, we have a well-motivated expectation for &l ; a = 
&L;f Ar/Sc, given that 8 C » 1.7, &L ; f ~ 1.4 (§ 2.1), and A r w 3 
from simulations (see § 3). Now, if r were independent of 
scale, then we could only measure a degenerate combination 
of the unknown quantities. However, a precise measurement 
of the power spectrum ratio can separate out the constant 
and BAO terms, thus yielding B\ and B2 separately, which 
in turn yields b r and the ratio &l ; a/&ii. 

Although it is implicit in the equations, r(k) and tlss 
are also declining functions of time. However, even at low 
redshift r(k) contains a signature of the BAOs, since the 
BAOs are still imprinted more strongly in the baryon fluc- 
tuations than in those of the dark matter or the total matter. 
This clear signature offers a chance to detect this effect, even 
if the various bias factors that we have assumed to be con- 
stant actually vary slowly with k. A detection of the effect 
can be combined with an estimate of 6 n from comparing P n 
with the underlying matter power spectrum (e.g., as mea- 
sured with weak lensing on large scales). Extraction of the 
value of 6l ; a would yield a new quantity in galaxy forma- 
tion, a combination of the way in which the luminosity of 
a galaxy depends on the baryonic content of its host halo, 
and of how this baryonic content depends on the underlying 
difference between the baryon and total density fluctuations. 

Another possibility is to compare the power spectra 
of luminosity density (or flux-limited number density) be- 
tween two different samples. Their ratio should again have a 
form similar to equation (|29[1 . from which the constant and 
BAO term can be separately measured. It i s well known that 
galaxy bias depends on galaxy luminosity (jLahav 199f|), but 
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Figure 1. The fractional baryon deviation r(k) = (<5 b /<5 tot ) — 1 
as a function of k, at various redshifts (2 = 0, 0.5, 1, 3, and 6, 
from top to bottom). 



here the bias would be scale dependent in a way that de- 
pends On Z/min- 



3 QUANTITATIVE PREDICTIONS 

For our quantit ative results, we u se the CAMB linear per- 
turbation code llLewis et alj|200d). with the WMAP 5-year 
cosmological parameters ( Komatsull2009h , matching the sim- 
ulation that we compare with below. 

We show the dependence of r on both wavenumber and 
redshift in Figure [T] At a given redshift, r(k) approaches a 
constant value (rLSs) at k 0.5 h/Mpc. Using tlss (itself 
a function only of redshift) we can separate out the two 
variables k and 2 in their effect on r, as shown in Figure [2] 
The function [r(k) /rxss] — 1 is independent of redshift (i.e., 
the curves for five different redshifts overlap very precisely) , 
so the k dependence of r is determined by a single, fixed 
function of k. Thus, the redshift dependence of r is the same 
at all k, and it suffices to show the dependence of j"lss- 
Figure [5] shows that, as noted in the previous section, t"lss 
indeed varies approximately in proportion to 1/a, but in 
detail the variation with redshift is slightly slower than that. 

As noted in the previous section, we expect the non- 
linear evolution that takes place during halo formation to 
magnify the gas depletion effect compared to the linear the- 
ory calculation. We can test this effect using the hydr ody- 
namical simulation o f lNaoz. Yoshida. fc Barkanal (|2010h . Al- 
though superficially it appears that they studied a quite dif- 
ferent regime (low-mass halos forming at high redshift), their 
results should be applicable here. In the linear theory, the 
gas depletion factor tlss is constant all the way from nearly 
the BAO scale (k ~ 0.5 h/Mp c) down to just above the 
Jeans scale (fc £ 100 h/Mpc). iNaoz. Yoshida. fc Barkanal 
investigated the gas depletion in virialized halos from 
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Figure 2. Top panel: [r(fc)/rLSs] — 1 as a function of k, at the 
same redshifts as in Figure [T] (the curves all lie on top of each 
other). Bottom panel: The quantity lOOrLss versus 1 + 2 (solid 
curve), or equivalently, the value of rLss in units of percent. Also 
shown is the function — 0.31/a (dotted curve). For both panels, 
in practice we set tlss = r (k = lh/Mpc). 

below the Filtering mass (which is a time-averaged Jeans 
mass) up to a thousand times higher mass scale. Thus, the 
most massive halos in their simulation were well into the 
large-scale structure regime, where pressure is negligible, 
and the effect we are interested in (i.e., non-linear gas de- 
pletion on large scales) should operate. 

Figure [3] shows that the fractional gas deple- 
tion measured in virializ e d ha los in the simulation of 
iNaoz. Yoshida. &: Barkanal (|201Cf ) was much larger than the 
depletion tlss predicted for linear perturbations at the 
halo virialization redshift z v i r - The simulated results can 
be reasonably fi10 either by multiplying 7~lss by a factor 
of 3.2, or by adopting tlss from a higher redshift z [where 
(1 + z) = 3.5(1 + Zvir)]- Additional simulations are required 
to test whether these results can indeed be extrapolated to 
our regime of much more massive halos at low redshift, but 
these results suggest that the gas depletion in halos is am- 
plified by a factor £ 3 compared to the linear regime. 

As an example of typical numbers, we consider an ex- 
ample with Cmin = D m i n = 0, b n = 2, and &L ; t = 1. As noted 
in § 12.41 we expect in this case &l ; a ~ 2.6, and also b r — 0.5, 
so in the observational ratio of equation (|29l) . B\ — 1.5 and 
B2 = 0.9. Thus, the oscillations in the square-rooted ratio 
of the luminosity and number density power spectra are at 
the level of 0.4% at z — 1 (measured from the first peak, i.e. 
at the lowest k, to the following trough; the variation from 
k = to the first peak is roughly twice as large). This is a 

2 There is a hint of a different slope with redshift in the simula- 
tion results compared to the fits. However, this mainly depends 
on a single point (at z v ; r = 12) and needs to be checked with 
further simulations. 



Figure 3. The fractional gas depletion in halos ver- 
sus redshift. We show the r esults from the simulation of 
INaoz. Yoshida. Barkana l|20ld ) (data points), where the cor- 
responding redshift (at which the virialized halos are identified 
in the simulation) is denoted z v ir- We compare the depiction as 
measured in the simulation to tlss a t ^vir (dotted curve), 3.2rLss 
at z v i v (long-dashed curve), and 7"lss at a value of (1 + z) equal 
to 3.5(1 + z v ir) (short-dashed curve). 

weaker effect by about a factor of five compared to the nor- 
mal BAOs in the total matter power spectrum. Thus, if high 
precision is achieved in the regular BAO measurement, then 
the scale-dependent bias that we have highlighted should 
also be measurable. 

This scale-dependent bias is unlikely to significantly af- 
fect the standard BAO measurements. Such measurements 
are usually carried out on the power spectrum of the galaxy 
number density. Scale-dependent bias enters this quantity 
only in proportion to C m i n (see equation I26[) . so it would 
be present only in a sample for which the flux limit plays 
a significant role. Even then, the effect on the BAO peak 
positions would be quite weak, since the BAOs in dtot are 
physically a result of the influence of the baryons on the 
dark matter. Thus, the peak positions in 5 to t and in Sb are 
nearly identical. For instance, even in the case that in equa- 
tion (126[l the coefficients (6 n + C* m in6L;t) and C m in6L;A are 
equal, the BAO peak positions are shifted by only ~ 0.3%. 



4 CONCLUSIONS 

We have shown that the variation in the baryon to mat- 
ter ratio imprinted by acoustic waves prior to cosmological 
recombination should result today in an oscillatory, scale- 
dependent bias of galaxies relative to the underlying matter 
distribution (see Figs. 1 & 2). The percent-level amplitude 
of this signature depends on how the typical luminosity of 
galaxies scales with the baryon mass fraction in the large- 
scale region in which they reside. Simulations suggest that 
this signature is significantly amplified by non-linear effects 
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during halo collapse (Fig. 3). The resulting amplitude can 
be measured from the ratio between the power spectra of 
fluctuations in the luminosity density and number density 
of galaxies (equation I29|l . An observational calibration of 
this amplitude would offer a new cosmological probe of the 
physics of galaxy formation. 

This effect may be marginally observable with current 
data, but it should certainly be feasible using future galaxy 
surveys (such as BOS^I or BigBOSSH). However, since the 
baryonic and the matter fluctuations have nearly identical 
BAO peak positions, the scale-dependent bias is unlikely to 
significantly affect the standard BAO measurements, even 
at percent-level precision. 
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